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Abstract 

We propose a general kinetic and hydrodynamic description of self-gravitating Brown- 
ian particles in d dimensions. We go beyond usual approximations by considering inertial 
effects and finite N effects while previous works use a mean-field approximation valid in a 
proper thermodynamic limit (N — * +oo) and consider an overdamped regime (£ — > +oo). 
We recover known models in some particular cases of our general description. We derive 
the expression of the Virial theorem for self-gravitating Brownian particles and study the 
linear dynamical stability of isolated clusters of particles and uniform systems by using 
technics introduced in astrophysics. We investigate the influence of the equation of state, 
of the dimension of space and of the friction coefficient on the dynamical stability of 
the system. We obtain the exact expression of the critical temperature T c for a multi- 
components self- gravitating Brownian gas in d = 2. We also consider the limit of weak 
frictions £ — > and derive the orbit-averaged-Kramers equation. 

1 Introduction 

Self-gravitating systems such as globular clusters and galaxies can be considered as a collection 
of N stars in gravitational interaction whose dynamics is described by the Hamilton equations 
of motion £Q. In statistical mechanics, this situation is associated with the microcanonical 
ensemble where the energy and the particle number are fixed [21 Ej- In a recent series of papers 
jH El El HI El El GDI EE] , we have proposed to consider a system of self-gravitating Brownian 
particles which are subject, in addition to the gravitational force, to a friction and a noise. Their 
dynamics is described by iV-coupled stochastic Langevin equations. In statistical mechanics, 
this situation is associated with the canonical ensemble where the temperature and the particle 
number are fixed • In previous papers, we have considered a mean-field approximation valid 
in a proper thermodynamic limit with N — > +oo and, in order to simplify the problem, we 
have studied a limit of strong friction £ — > +oo, or equivalently a large time regime t ^> 
In these approximations, the problem is reduced to the study of the Smoluchowski-Poisson 
(SP) system. We have also introduced a generalized class of stochastic processes and kinetic 
equations in which the diffusion coefficient (or the friction/mobility) is allowed to depend on the 



concentration of particles fHlEJEHElE]- This can model microscopic constraints (e.g. close 
packing) acting on the particles when their local concentration becomes large. The evolution of 
the system is then described by the generalized Smoluchowski-Poisson (GSP) system involving 
a barotropic equation of state p(p) specified by the stochastic model. This mean-field nonlinear 
Fokker-Planck (MFNFP) equation admits a Lyapunov functional, determined by the equation of 
state, which can be interpreted as a generalized free energy [T3]. Thus, this model is associated 
with a notion of (effective) "generalized thermodynamics" in /i-space. In the classical case 
where the diffusion coefficient is constant, we recover an isothermal equation of state p = 
phsT/m associated with the Boltzmann free energy but more general equations of state can 
be considered. Interestingly, the same type of drift-diffusion equations are encountered in 
mathematical biology to describe the chemotactic aggregration of bacterial populations, in 
relation with the Keller-Segel model ^SIEHH^O]- The analogy between self-gravitating Brownian 
particles and bacterial populations is developed in p. 

Here, we propose to generalize these models so as to take into account finite N effects 
and inertial effects (finite friction £). We thus propose a general kinetic and hydrodynamic 
description of self-gravitating Brownian particles starting directly from a system of N coupled 
Langevin equations of motion with long-range attractive interactions. We shall extend the 
technics developed in astrophysics to our problem of self-gravitating Brownian particles. In 
particular, we shall derive the appropriate expression of the Virial theorem for these systems 
and study their linear dynamical stability by a method similar to that developed by Eddington 
[2*T] and Ledoux |22] f° r barotropic stars described by the Euler equations. We shall make the 
link between parabolic and hyperbolic models by considering an intermediate model taking into 
account the inertia of the particles as well as a friction force. The Euler equations are recovered 
for £ = and the Smoluchowski equation is obtained in the limit £ — > +00. 

The paper is organized as follows. In Sec. |2l we introduce general kinetic and hydrodynamic 
models of self-gravitating Brownian particles starting from coupled Langevin equations. We 
derive the iV-body Fokker-Planck equation ( subsection 12 .lj) . the mean- field Kramers equation 
( subsection I2.2JI and the generalized mean-field Kramers equation ( subsection 12. 3|) . Then, we 
take the hydrodynamic moments of these equations and derive the damped Jeans equations 
(subsection I2.4J1 and the damped barotropic Euler equations (subsection I2.5J1 by closing the hi- 
erarchy of moments with a local thermodynamical equilibrium (L.T.E.) hypothesis. We obtain 
the mean- field Smoluchowski equation in a strong friction limit £ — > +00 and, in subsection 12. bl 
we derive the orbit-averaged-Kramers equation is a weak friction limit £ — > 0. In Sec. EJ we es- 
tablish the general expression of the Virial theorem for self-gravitating Brownian particles from 
the damped Jeans equations (subsection 13. lj) and from the damped Euler equations (subsec- 
tion (sT2J). We also consider the effect of correlations due to finite N effects ( subsection 13 .3|) . In 
Sec. El we study the linear dynamical stability of an inhomogeneous stationary solution of the 
damped barotropic Euler-Poisson system and investigate the effect of the friction coefficient on 
the evolution of the perturbation. In Appendix [X] we give a short complement concerning the 
stability of polytropic systems and in Appendix El we derive the exact expression of the Virial 
theorem starting directly from the stochastic equations of motion. We show that the Virial 
theorem takes a very simple form in dimensions d = 2 and d — 4 and analyze the consequences 
of this simplification. In Appendix O we study the linear dynamical stability of homogeneous 
stationary solutions of the damped barotropic Euler equations (for an arbitrary potential of 
interaction) and obtain a generalization of the Jeans instability criterion. In the Conclusion, 
we discuss the different regimes in the evolution of Hamiltonian and Brownian systems with 
long-range interactions ^2] distinguishing the phase of violent collisionless relaxation, the col- 
lisional evolution due to finite N effects and the "collisional" evolution due to imposed friction 
and stochastic forces for Brownian systems. 
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2 Kinetic and hydrodynamic models of self-gravitating 
Brownian particles 

The Smoluchowski-Poisson system that has been extensively studied in previous papers jU 
El El 13 El El [Jill UJ] describes a gas of self-gravitating Brownian particles in a mean-field 
approximation (valid for iV — > +00) and in a strong friction limit £ — > +00. In this section, we 
introduce more general kinetic and hydrodynamic models of self-gravitating Brownian particles 
that go beyond these approximations. 



2.1 The iV-body Fokker- Planck equation 

Basically, a system of self-gravitating Brownian particles is described by the iV coupled stochas- 
tic equations 



—jjT = -£vj - raVi£/(ri, 



T N ) + y/2DRi(t), 



(2) 



where — £vj is a friction force and Rj(£) is a Gaussian white noise satisfying (Rj(t)) = 
and (R at i(t)Rb t j(t')) = 5ij5 a bS(t — t'), where a,b = l,...,d refer to the coordinates of space 
and i,j = 1,...,N to the particles. We define the inverse temperature (3 = l/k B T through 
the Einstein relation £ = D/3m (see below). For £ = and D = 0, Eqs. (UJ)-© reduce to 
the Newton-Hamilton equations of motion describing the ordinary self-gravitating gas with a 
Hamiltonian 



N 1 



(3) 



where U(ri, ...,r^v) = ~ r i) an< ^ M ( ri — r i) = — G/[(d — 2)|r^ — r^- 1 2 )] denotes the 

gravitational potential of interaction in d dimensions (w(rj — r^) = Gin |r$ — | for d = 2). In 
this paper, we shall be particularly interested by the gravitational interaction, but we stress 
that our formalism remains valid for a more general class of binary potentials of interaction of 
the form u(\t{ — rj\). The evolution of the iV-body distribution function is governed by the 
Fokker-Planck equation 12J : 



dP, 



N 



dt 



N 

E 

i=l 



dP, 



N 



Or, 



dP, 



N 



dvi 



N 



_d_ 



(4) 



In the absence of friction and diffusion (£ = D = 0), it reduces to the Liouville equation. 
The Liouville equation conserves the energy (E) = J P^H FJ^ dr^dwi and the Gibbs entropy 
S = —k B f P N In P N dYidv\. . .dr jydv N (more generally, any functional of P^) defined on T-space. 
This corresponds to a microcanonical description. Alternatively, in the Brownian model, the 
temperature T is fixed (instead of the energy) and the Fokker-Planck equation (J1J) decreases 
the Gibbs free energy F = (E) — TS which can be written explicitly, 



F[P N ] = J PNHYldridVi + kgT J PtflnPjvJJdridvt. 



(5) 
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This corresponds to a canonical description. One has 



/ 1 / dP \ 
F=-J2J Jp^ ( D ~d^r + £ P ^J dr 1 dv 1 ...dr N dv N < 0. (6) 

Therefore, the free energy plays the role of a Lyapunov functional for the iV-body Fokker- 
Planck equation. At equilibrium, F = implying that the r.h.s. of Eq. (0J vanishes. The l.h.s. 
(advective term) must also vanish independently. From these two requirements we find that 
the stationary solution of Eq. (J3J is the canonical distribution 1 

P N (n, vi, r N , y N ) = J-e-^Ef^ i+mu^,...^)) ^ (7) 

provided that the coefficients of diffusion and friction are connected by the Einstein rela- 
tion £ = Df3m. The partition function Z(f3) is determined by the normalization condition 
J P/v n« dvid~Vi = 1. The canonical distribution (JJJ) minimizes the free energy F at fixed parti- 
cle number. Introducing the reduced probability distributions 



F,(xi, ...,Xj) = J P N (x 1 ,...,x N )dx 1i+ i...dx N , 



where x = (r, v), we can readily write down a hierarchy of equations for Pi, F 2 etc. The first 
equation of the hierarchy is 



at oyi J avi avi 



dvi 



(9) 



where F(2 — > 1) = —mduu/dri = Gm{v2 — ri)/|r2 — ri| d is the force (by unit of mass) created 
by particle 2 on particle 1. Note that this equation is exact (i.e. valid for all N) and takes 
into account the correlations between the particles encapsulated in the two-body distribution 
function P 2 . For D = £ = 0, we recover the first equation of the BBGKY hierarchy. We shall 
give the form of the Virial theorem associated to Eq. Q in Sec. 13.31 Before that, we consider 
the mean-field limit of this equation valid for N — > +oo. 



2.2 The mean-field Kramers equation 



In a properly defined thermodynamic limit [12] , we can show that the cumulant of the two-body 
correlation function is of order 1/N. Therefore, for iV — > +oo, we can implement the mean-field 
approximation 

P 2 (n, vi, r 2 , v 2 , t) = Pi(ri,vi, t)P x {r 2 , v 2 , t) + 0(1/N). (10) 

Substituting this result in Eq. © and introducing the distribution function / = NmP x , we 
obtain the mean-field Kramers equation 



df df ^ df d 
+ v • — — V$ ■ — = — 
dt dr dv dv 



av 



order to properly define a strict statistical equilibrium state for self-gravitating systems, one has to 
introduce a small-scale regularization otherwise F[Pn] has no minimum and Eq. Q is not normalizable. Thus, 
in Eq. (JJJ, it is implicitly understood that U is a regularized potential. Note that physical statistical equilibrium 
states unaffected by the small-scale regularization exist in the form of long-lived metastable structures that are 
local minima of the Boltzmann mean- field free energy Fb [f] defined in Eq. I) bill . We refer to [21] for a physical 
discussion of these issues. 
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where <3>(r, t) — J u(r — r')p(r', t)dr'. This equation is non-local because the potential $(r, t) is 
induced by the density p(r, t) = J fdv of the particles composing the whole system (it is not an 
external potential). Thus, for self-gravitating systems, Eq. (|11|) has to be solved in conjunction 
with the Poisson equation 

A$ = S d Gp. (12) 

In the absence of friction and diffusion (D = £ = 0), the mean- field Kramers equation reduces to 
the Vlasov equation lj. The Vlasov equation conserves the energy E = | J fv 2 drdv + ^ J p§dr 
and the Boltzmann entropy S B = — f(f /m)hi(f /m)drdv (more generally all the functionals 
of / called the Casimirs) defined on p-space. Alternatively, the Kramers-Poisson (KP) system 
(|11 )) -([12 )) involves a fixed temperature and decreases the Boltzmann free energy F B = E — TS B 
which can be written explicitly 



F B [f] = - [ fv 2 dvdv+- [ p§dr + k B T [ — \n—drdv. 
2 J 2 J J m m 



(13) 



Indeed, one has 



B 



h{ D ^ + ^ fv ) drdv -°- (14) 



At equilibrium, F B = implying that the r.h.s. of Eq. Ijlljl vanishes. The l.h.s. (advective 
term) must also vanish independently. From these two requirements, and using the Einstein 
relation, we find that the stationary solutions of the Kramers-Poisson system fjllj) correspond 
to the mean-field Maxwell-Boltzmann distribution 

/ = Ae -Pm{4+Hr)) ( 15 ) 

which has to be solved in conjunction with the Poisson equation ()12p. The stable mean-field 
Maxwell-Boltzmann distribution minimizes the Boltzmann free energy F B [f] at fixed mass. It 
is both thermodynamically stable (in the canonical ensemble) and linearly dynamically stable 
with respect to the KP system ^3] ■ We note that the equilibrium one-body distribution function 
(J 15)) can be obtained from the iV-body canonical distribution (J7J) by constructing an equilibrium 
BBGKY-like hierarchy and implementing a mean-field approximation [T2] . On the other hand, 
the Boltzmann free energy (|13|) can be deduced from the free energy of the iV-body system (J3J) 
by making the mean-field approximation P/v(l, -..,N) = Yli^i (0 

2.3 The generalized non-local Kramers equation 

For sake of generality, we shall consider the case where the diffusion coefficient explicitly depends 
on the distribution function. Thus, in Eq. ([lip, we set D(f) = DfC"(f) where C is a convex 
function, i.e. C" > 0, and D is a constant. In that case, we obtain the generalized mean-field 
Kramers equation [THj: 



df df ^ df d 

— + v • — - V$ • — = 

dt dr dv dv 



DfC"(f)^+Zfv 



(16) 



This equation can be obtained in the mean-field limit of a generalized iV-body Fokker-Planck 
equation associated with a generalized class of stochastic processes of the form ([T|)-(j2"]) where 
the diffusion coefficient depends on /(rj,Vj,i) For C = /In/, we recover the classical 
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Kramers equation pip. However, Eq. p(j|) can describe more general situations such as quantum 
particles with exclusion or inclusion principles (fermions, bosons, quons), lattice models, non- 
ideal effects etc... These generalized Fokker-Planck equations are associated with an effective 
thermodynamical formalism (ETF) in //-space [13]. In particular, the generalized Kramers- 
Poisson (GKP) system decreases the free energy 

F\f]=E-TS= f fjdrdv+~f p^dr + T Jc(f)drdv, (17) 

where the last term can be interpreted as a generalized entropy S = — J C(f)drdv and we have 
defined the effective temperature T = 1/(3 through the relation £ = D/T (effective Einstein 
relation). One has 

F = ~ S Tz( D f C V)7r + ZfA d ^ < °- (18) 



At equilibrium, we find that the stationary solutions of the generalized Kramers equation (JTSj) 
are determined by the integro-differential equation 

C(J) = -p(j + *\ -a, (19) 

where $(r, t) = J u(r — r')/(r', v', t) dr'dv'. Since C is convex, this relation can be inversed 
to give / = F(/3e + a) where F(x) = (C") _1 (-a;) and e = ^ + $ ( r )- We note tnat tne 
equilibrium distribution determined by Eq. (fT9~)) is a function f — f(e) of the individual energy 
e alone which is monotonically decreasing (for (3 > 0). This equilibrium distribution function 
extremizes the free energy (fT7|) at fixed mass. Furthermore, only a (local) minimum of free 
energy is linearly dynamically stable with respect to the GKP system |13j . 

2.4 The damped Jeans equations 

We shall now determine the hierarchy of moment equations associated with the generalized 
Kramers-Poisson system. Defining the density and the local velocity by 

P = J fdv, pu = J fvdv, (20) 

and integrating Eq. (fTtjj) on velocity, we get the continuity equation 

f t + V • (pu) = 0. (21) 

Next, multiplying Eq. p6p by v and integrating on velocity, we obtain 

^(piH) + Wj {pWj) = -^7 " Pq^ ~ ^ ( 22 ) 
where we have defined the pressure tensor 

Pij = / fwiWj dv, (23) 



where w = v — u is the relative velocity. In the absence of diffusion and friction (D = £ = 0), 
we recover the Jeans equations of astrophysics that are derived from the Vlasov equation [T]. 
For self-gravitating Brownian particles, the equivalent of the Jeans equations (|2ip - (|22p include 
an additional friction force — £u. Using the continuity equation, Eq. f)22|) can be rewritten 



[dm dui\ dPa 9$ „ 

p [w + ^) = --d^- p te i -t pUi - (24) 
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2.5 The damped barotropic Euler equations 

By taking the successive moments of the velocity, we can obtain a hierarchy of hydrodynamical 
equations. Each equation of the hierarchy involves the moment of next order. The ordinary 
Jeans equations that are based on the Vlasov equation are difficult to close because the Vlasov 
equation admits an infinite number of stationary solutions. Therefore, a notion of thermody- 
namical equilibrium is difficult to justify in the usual point of view (see, however, in the 
context of the theory of violent relaxation). In the present case, the situation is simpler because 
the Kramers equation admits a Lyapunov functional (fTTj) and a unique stationary distribution 
defined by Eq. (|19|) . If we are sufficiently close to equilibrium, it makes sense to close the 
hierarchy of equations by using a condition of local thermodynamic equilibrium. We shall thus 
determine the pressure tensor Eq. (j2HJ) with the distribution function /l.t.e defined by the 
relation 

,2 



C'(h.T.E.) = -P 



W 

T + A(r,t) 



(25) 



This distribution function minimizes the generalized free energy Eq. (fTTj) at fixed temperature 
T, local density p(r, t) and local velocity u(r, t). The function A(r, t) is the Lagrange multiplier 
associated with the density field. It is determined by the requirement 

p(r,t) = J f L .T.E.dv = p{\{r,t)]. (26) 

At equilibrium, we recover the distribution function Eq. (|19|) with u(r) = and A(r) = $(r) + 
a/f3. Using the condition Eq. (|25|) of local thermodynamic equilibrium (LTE), the pressure 
tensor Eq. (}2*3*j) can be written = p5ij with 



P(r,t) = ~ J Jl.t.e.w 2 dw = p[\(r,t)]. 



(27) 



The pressure is a function p = p{p) of the density which is entirely specified by the function 
C(f), by eliminating A from the relations Eq. ()26|) and Eq. (|27|). We note furthermore that, 
using Eq. (j23|) and integrations by parts, the previous equations (f2lj|) and (|2T|) easily lead to 
= \ J %w 2 dw = ±VA / f£wdw = -VA / fdw = -pVA, hence p'{p) = -pA'(p). In the 
case of Brownian particles described by the ordinary Kramers equation (Jll)) with C(f) = f In /, 
the equation of state determined by Eqs. ()26|) and (}2~Tj) is the isothermal one p = ^-p- More 
generally, we obtain the damped Euler equations for a barotropic gas [T3] : 

^ + V • {pa) = 0, (28) 

^ + (u-V)u = — Vp-V$-£u. (29) 
at p 

These equations decrease the free energy 

F[p, u] = J p J'^l dp'dv + p$dr + / P y dr, (30) 

which can be deduced from the free energy (fTTjl by using the local thermodynamic equilibrium 
condition (|23|) to express F[f] as as a functional of p and u, using F[p, u] = ^[/l.t.b.] (see [Ej 
for details). A direct calculation yields 

F = -£ pu 2 dr < 0. (31) 
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At equilibrium, F = implying u = 0. Then, Eq. (|29|) yields the condition of hydrostatic 
balance 



Vp + pV$ = 0, (32) 

which also results from Eq. (|19|) . Indeed, for / = /(e) with e = v 2 /2 + $(r), one has p = 
f /(e) dv, P = \j f{t)v 2 dv so that 

Vp = i / /'(e) V$v 2 dv = iv$ j y~-vdv= -V$ y /dv = -pV$. (33) 

The damped barotropic Euler equations (j2E|) - (|2n|) are interesting as they connect hyperbolic 
models to parabolic models. Indeed, for £ = we recover the standard barotropic Euler 
equations describing the dynamics of gaseous stars [H 123 |2H 121] or the formation of large-scale 
structures in Cosmology j2Hj. Alternatively, in the strong friction limit £ — > +oo, we can neglect 
the inertial term in Eq. (J2"H|) and we obtain 

u = -^-(Vp + pV$) + 0(£" 2 ). (34) 

Substituting this drift term in the continuity equation (|28|). we get the generalized Smoluchowski 
equation 



This equation decreases the free energy 



i(Vp + pV$) 



(35) 



F[ P ] = fpfyr d P' dr + \j (36) 

which is obtained from Eq. (|3U|) by neglecting the last term of order 0(£ -2 ). A direct calculation 
yields 

— (Vp + pV$) 2 dr < 0. (37) 

£p 

It should be recalled that the damped Euler equations (|2"%|) - (f2"H|) remain heuristic because 
their derivation is based on a Local Thermodynamic Equilibrium (L.T.E.) assumption (j2Sj) 
which is not rigorously justified. However, using a Chapman- Enskog expansion, it is shown 
in J7| that the generalized Smoluchowski equation (J33j) is exact in the limit £ — > +oo (or, 
equivalently, for times t ^> £ _1 ). The generalized Smoluchowski equation can also be obtained 
from the moments equations of the generalized Kramers equation by closing the hierarchy, using 
£ — ► +oo (see Sec. 9 of [To]). 



2.6 The orbit-averaged-Kramers equation 

We shall consider here the opposite limit of low friction £ — > 0. In the case where the term 
in the r.h.s of Eq. (fTT|) is small, we can obtain a simplified equation for the evolution of the 
distribution function by averaging the kinetic equation over the orbits of the particles. In 
the case of Hamiltonian self-gravitating systems described by the Landau equation, this leads 
to the orbit-averaged-Landau equation introduced by Henon j2H]- We shall here derive the 
orbit-averaged-Kramers equation for self-gravitating Brownian systems. 
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Let us first rewrite the mean-field Kramers equation (Jll|) in the form 



| +v .|_ V<t .|. W , ( 3S) 

We consider the case where £ — > for fixed (3 so that £ ~ D ~ and the operator can 
be considered as small. If we take Q(f) = 0, we obtain the Vlasov equation. We shall assume 
that the system has reached a stable stationary distribution of the Vlasov equation of the form 
f — /(e) which depends only on the energy e = \ + <3>(r, t) of the particles. This is a particular 
case of the Jeans theorem for spherical systems. Such steady solution can arise from a process 
of violent collisionless relaxation [3U1[T2]. Since Q(f) ^ 0, the distribution function will change 
due to the terms of friction and diffusion that are present in the stochastic equation Eq. © 2 . 
However, if £ — *• 0, this change will be slow so that the latter forces cause only a small variation 
on the energy. We shall therefore consider that f(r, v,t) ~ /(e, t) remains a function of the 
energy alone that slowly evolves due to imposed stochastic forces. Noting that 

d df d§df 

M r,,, M = -t + £ (39) 



we can rewrite the kinetic equation (JHHj) in the form 

d -l + d A d _L = Q (f) (40) 

dt dt de WJ ' 1 ' 

Since / depends only on the energy, the system is spherically symmetric. Then, the phase space 
hypersurface with energy smaller than e is 

q(e,t) = 16tt 2 / r 2 drv 2 dv = -— / [2(e - $)] 3/ Vdr, (41) 

where r max (e, t) is the largest radial extent of an orbit with energy e at time t. It is determined 
by the equation e = <&{r maxi i) corresponding to v = 0. The previous relation can be written 
more compactly as 

q(e,t) = — / v 3 r 2 dr, (42) 



where v = a/2(c — $(r, t)). The phase space hypersurface g(e,t)de with energy between e and 
e + de is given by 

g (e, t) = ^- = 16vr 2 / [2(e - <$>)} l ' 2 r 2 dr = 16tt 2 / vr 2 dr. (43) 

d e Jo Jo 

Now, the density of particles in the hypersurface between e and e + de is uniform since the dis- 
tribution function depends only on the energy. In fact, the system evolves on a short timescale 
~ to by purely inertial effects (corresponding to the advective terms in the l.h.s of Eq. (jHEI)) so 

2 In the case of Hamiltonian systems where £ = D = 0, the distribution function changes due to finite N effects 
representing close encounters between stars. These encounters are usually modelled by the Landau operator |12| 
which is of order 1/N <C 1. On the other hand, for self-gravitating Brownian particles the evolution is modeled 
by the Kramers operator (obtained for N — * +oo) which is of the order (to (where to is the dynamical time). 
The comparison between the evolution of Hamiltonian and Brownian systems with long-range interactions is 
further discussed in the Conclusion of this paper. 



9 



as to establish this quasi-stationary regime where / ~ f(e,t). We shall therefore average the 
kinetic equation (|4*0*)l on each hypersurface of iso-energy using 



jj-maz Xvr 2 dr 
(X)(e,t) = rrmax — — — 
J vr z dr 

for any function X(r, v,t). Thus, the orbit-averaged kinetic equation can be written 

'df d§ df 



(44) 



/ ' max 

16n 2 / r 2 dr v 
Jo 



3t + mt- Q(!) 



0. 



The first term in bracket is 

167T 

Since 



max 

2 / „2 



r dr v- = g(e,t)- = --. 



dq 1R 2 f rmax <9$ 2 



.2 / ^2 



the second term in bracket can be written 

16tt 

Finally, since 



; ., (/r ^ = _dqd£ 
dt de dt de ' 



vQ(f) = D§- e 



.3 {df 



the last term in bracket is 
16tt 



max 

2 / „2 







(9 

r a dr uQ(/) = 3D— 



Regrouping all these results, we obtain the orbit-averaged-Kramers equation 



dqd£ _ dqd£ = 3£) _9 
de dt dt de de 



' ~ie + 



(45) 



(46) 



(47) 



(4£ 



(49) 



(50) 



(51) 



167T 2 r r max(t,t) 

q( e ,t) = —j [2(e-$(r,t))] 3 /V 2 rfr, 



(52) 



r — ) = 16vr G 



f(e,t)[2(e-$(r,t))] 1/2 r 2 de, 



*(»•,*) 



(53) 



where the last equation is the Poisson equation. It is easy to verify that the free energy is 
monotonically decreasing (F < 0) and that the stationary solution of Eq. (|51|) is the Boltzmann 
distribution / = Ae~^ me . This equations will be studied in a future communication. Note also 
that in d = 1, q = § vdx, g = 2 § dx/v and in d = 2, q = 2tt 2 J Q rm v 2 rdr and g = 2ii 2 r 2 n . 
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3 The Virial theorem for Brownian particles 



3.1 The Virial theorem from the damped Jeans equations 

We shall give here the form of the Virial theorem appropriate to the damped Jeans equations 
fj21|l -(|22 )l . The only difference with the standard Jeans equations is the presence of a friction 
term. We shall thus only give the final result and refer to pQ for the details of calculation. The 
damped Virial theorem can be written 

ll& + \^ = 2KlJ + Wij ~ \ f {PikXj + PjkXi) dSk - (54) 

We have included boundary terms which must be taken into account if the system is confined 
within a box. The tensor of inertia and the potential energy tensor W%j are defined in Paper 
I. The kinetic energy tensor is defined by 

K ti = n / f v i v 3 dw - ( 55 ) 



It can be written as 



where 



2 



K a = T ij + ( 56 ) 



Tij = ^ J pUiUj dr, Uij = J Pij dr. (57) 

Note that the tensors and depend on the distribution function /(r, v, £), not only on 
hydrodynamic variables. The scalar Virial theorem takes the form 

~ + =2K + W«-f P ik x % dS k , (58) 

where / is the moment of inertia and 

K = ^Jfv 2 dvdv, (59) 

is the kinetic energy. It can be written 

K = T + ±I1, (60) 

where 

T = ^Jpu 2 dr, U = djpdr. (61) 

In the absence of diffusion and friction (D = £ = 0), we recover the usual expression of the 
Virial theorem issued from the Jeans equations pQ. For Brownian particles, the novelty is the 
presence of a damping term £J. 

As pointed out in Paper I, the moment of inertia depends on the origin O of the system of 
coordinates. Let R(i) = {1/M) J prdr denote the position of the center of mass with respect 
to the origin O. Using the equation of continuity (|28j) . we find that MdH/dt = P where 
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P = J pudr is the total impulse. Using the Jeans equation (J22|) . we find that dP/dt = — £P. In 
our case, there exists an absolute referential {7V). Indeed, in writing Eq. (0) we have implicitly 
assumed that our Brownian particles evolve in a fluid that it motionless. Otherwise, the friction 
force in Eq. (j2J) has to be modified according to — £(Vj — U) where U is the velocity of the 
fluid |T7j. We must work therefore in this referential {7V}- If we denote by Po the initial value 
of the impulse, we get P(t) = Poe - ^. If now Ro denotes the initial position of the center of 
mass with respect to O, we find that 

R(t) = R + i|(l-e^). (62) 

Therefore, at t — > +00 the center of mass has been shifted by a quantity P /M£. In the strong 
friction limit £ — > +00, we find that the center of mass is motionless (Paper I). 



3.2 The Virial theorem from the damped Euler equations 

The Virial theorem associated to the damped barotropic Euler equations (|28| ) - (j29|) can be 
deduced from Eq. (J53J) by using the fact that = p(p)5ij. This yields 

\ d ~di? + h^ikt = 2Tij + \ mii + Wij ~\j p ^ ikXj + 6jkXi ^ dSk ' (63) 



where each quantity is now expressed in terms of hydrodynamic variables. At equilibrium, if 
no macroscopic motion is present (T = 0) and if we can neglect the boundary term, we get 

Wy = --n<%. (65) 

In the absence of diffusion and friction (D = £ = 0), we recover the usual form of Virial 
theorem issued from the barotropic Euler equations 1J. Alternatively, in the strong friction 
limit £ — > +00, we can neglect the term I in front of £/. Furthermore, since the velocity field 
scales as u = 0(£ -1 ), the kinetic energy tensor Tjj = 0(£~ 2 ) can also be neglected. Therefore, 
the Virial theorem associated with the generalized Smoluchowski-Poisson system (jHB]) - (fT2*jl can 
be written 

1 dl 1 f 



2 dt 

where 



1 ^ = 2K+W li - <b pr-dS, (67) 



K i:j = ^KSij, K = ^ / pdv. 



is the expression of the kinetic energy to leading order in the limit £ — > +00 where u(r, t) can 
be neglected. We thus recover the results of Paper I starting directly from the GSP system. 
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3.3 The effect of correlations 



If we account for the effect of correlations (due to finite values of N) between the particles and 
use the exact kinetic equation Q, we obtain the exact damped Jeans equations 

d -l + V ■ (pu) = 0, (69) 



Q $ Q P- ■ f 'x! ■ x ■ 

-^Oi) + far(P u * u ^ = —fa; + Gm2 J | r /_ r |> 2 ( r ' r '>*) dr ' ~ ^> ( 7 °) 



dt dxj dxj 

where we have introduced the spatial correlation function 



p 2 fa,r 2 ,t) = N(N - 1) / P 2 (ri, vi, r 2 , v 2 , t) dv x dv 2 . 



(71) 



In the mean-field approximation p 2 (l,2) = p(l)p(2), we recover the damped Jeans equations 
(I21j) - (|22j) . The Virial theorem is now given by 

Yw~ + Y~lt = 2K ^ + W ^ r ~ \ f {PlkX > + Plkxi] dSk ' (72) 

where 

Urcorr Gm 2 f fa ~ x\) fa ~ x'j) 

w corr = __ J p2 (r jT ) | r _ r ,| d drdT > ( 73 ) 

is a generalization of the potential energy tensor accounting for correlations between particles. 
In the strong friction limit £ — > +oo, the Virial theorem reduces to 

\^ = k J pdr + W™ rr - l - ^p(5 ikXj + djtxi) dS k . (74) 

If we consider the case of Brownian particles with an isothermal equation of state p = pk B T /m 
and if we focus on a space with d = 2 dimensions where 

W corr _Gm 2 f , „ , „,„ ^Gm 2 

u 2 

the scalar Virial theorem takes the form 



I" Gm 2 
J P2 (y,v')dvdv' = -N{N -l)^, (75) 



with a critical temperature 



l^ = 2Nk B (T-T c )-2PV, (76) 



Gm 2 {N-l) 

k BT c = . (77) 



These results are valid for an arbitrary number of particles. For N ^> 1, using iV — 1 ~ N, we 
recover the results of Paper I. 
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4 Dynamical stability of self-gravitating Brownian par- 
ticles 



We shall now investigate the linear dynamical stability of a stationary solution of the damped 
barotropic Euler-Poisson system (f2"g )) -(|2"9 j) satisfying the condition of hydrostatic balance (|32|). 
We shall determine in particular the equation of pulsations satisfied by a small perturbation 
around this equilibrium state. We shall investigate the influence of the friction parameter £ 
on the pulsation period and make the connection between the standard Euler-Poisson system 
(hyperbolic) obtained for £ = and the generalized Smoluchowski-Poisson system (parabolic) 
obtained for £ — > +00. The linearized damped Euler-Poisson equations are 

^ + V ■ (p5u) = 0, (78) 

p-^ = -V(p'(p)Sp) - pV5$ - <5pV$ - £p5u, (79) 

A5$ = S d G5p. (80) 

Considering spherically symmetric systems and writing the evolution of the perturbation as 

5p ~ e xt , we get 

Atfp+-i r J : (r*-V«) = 0, (81) 

X P 5u = -%'{ P )5p) -P^-- 6pf- - £p5u, (82) 
dr dr dr 

1 d =S&p. (83) 



r d 1 dr\ dr 
As in Paper I, we introduce the function q(r) defined by 



" = (84) 



The continuity equation then yields 



After some elementary transformations similar to those of Paper I, Eq. (|82|) can be put in the 
form 

d_ ( y\p) dq\ + Gq_ = A(A + Q 
dr\Sdpr d ~ 1 dr] r d ~~ l Sdpr^ 1 

The case of barotropic stars described by the Euler-Poisson system corresponds to £ = 
PU 123 1213 I2Z|- The case of self-gravitating Brownian particles described by the generalized 
Smoluchowski-Poisson system is recovered for £ ^> A (see ^3] and Paper I). We can therefore 
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use the results of Paper I by making the substitution £A — ► A(A+£). Therefore, an approximate 
analytical expression for the eigenvalue A is given by 

A(A + = (cPy + 2-2d)(d-2)y, (d ^ 2) (87) 

A(A + = -(7-l)^-, (d=2). (88) 

The friction coefficient £ affects the evolution of the instability but it does not change the 
instability threshold (determined by the sign of the l.h.s. of Eqs. (|87| ) -([88| ) ). The unstable case 
corresponds to A(A + £) = cr 2 > 0. The two eigenvalues are 

A ± = ± vf^ . (89) 

Since A + > 0, we see that the perturbation grows exponentially rapidly as e x+t . The stable 
case corresponds to A(A + £) = — cr 2 < 0. The two eigenvalues are 



A± = (90) 

If £ 2 — 4<T 2 > 0, then A± < and the perturbation decreases exponentially rapidly without 
oscillating. This is the case in particular for Brownian particles described by the Smoluchowski 
equation (£ — > +oo) for which A = — o" 2 /£ (Paper I). Alternatively, if £ 2 — 4<r 2 < 0, then 
A± = (— £ ± i\jAo 2 — £ 2 )/2 and we have slowly damped oscillations with a pulsation u> = 
|y / 4cr 2 — £ 2 and a damping rate £/2. This is the case in particular for barotropic stars (£ = 0) 
which oscillate with pulsation = cr without attenuation. The separation between these two 
regimes (pure damping vs damped oscillations) is obtained for £ = 2a at which u = 0. This 
suggests to introducing the dimensionless parameter 

F ^M&TY (91) 

measuring the efficiency of the friction force. The critical values are F = and F = —4. If 
F < —4 the system is stable and a perturbation is damped out exponentially rapidly without 
oscillating. If —4 < F < the system is stable and a perturbation exhibits damped oscillations. 
The pulsation vanishes for F = —4 and the damping rate vanishes for F — 0. For F > 0, the 
system is unstable. Using Eqs. (f8~7 ]) -([SS j) . the parameter defined in Eq. (|5T|) is explicitly given 
by 

1 i 2 I 



(d = 3) (92) 
(d = 2) (93) 



U 2 J 



(37 - 4) W 

-1 £ 2 / 
7 - 1 CM 2 ' 



Dimensionally, this parameter scales as |F| ~ £ 2 R d /GM. It can also be written \F\ ~ (£^) 2 
where ~ 1/ ^/pG is the dynamical time pQ . The dynamical stability of a homogeneous system 
(for a general form of potential of interaction) is treated in Appendix O 
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5 Conclusion 



In this paper, we have introduced general models of self-gravitating Brownian particles (stochas- 
tic iV body, kinetic, hydro dynamic,...) that relax the simplifying assumptions that are usually 
considered: mean-field approximation for N — > +00 (thermodynamic limit) and overdamped 
approximation for £ — > +00 (strong friction limit). These general models show the connection 
between previously considered models and offer a unifying framework to study these systems. 
We have focused here on the case of self-gravitating systems but most of our results also apply 
to the problem of chemotaxis in biology. This will be specifically considered in another paper 
where we discuss inertial models of bacterial populations. 

It should be emphasized that the Brownian model (fTj)-(|2|) contains the standard Hamiltonian 
model of stellar dynamics |1 as a special case since the Langevin equations reduce to the 
Hamilton equations for £ = D = 0. We expect therefore to have different regimes depending 
on the value of the parameters. To characterize these regimes properly, it is useful to introduce 
different timescales: (i) the dynamical time to ~ 1/y/pGl (Kepler time) is the typical period 
of an orbit or a typical free-fall time pQ (ii) the collisional relaxation time tu ~ (N/ In N)tu 
(Chandrasekhar time) is the typical time it takes a stellar system (Hamiltonian) to relax to the 
Boltzmann distribution e~ /3me due to close encounters. This relaxation is due to finite N effects 
PT] (iii) the friction time tg ~ (Kramers time) is the typical time it takes a Brownian 
system to thermalize, i.e. to have its velocity distribution close to the Maxwellian eT^ mv I 2 
[32] . This thermalization is due to the combined effect of imposed friction and diffusion in 
the Langevin model (fT|)-P|). It is due to a thermal bath (of non-gravitational origin) not to 
collisions (finite N effects). We can now distinguish different cases: 

(1) The case to <C t^ <C ts (£ — > 0) corresponds to Hamiltonian systems. For t <C £r, 
the system is described by the Vlasov-Poisson system. There is first a phase of violent colli- 
sionless relaxation on a timescale ~ to leading to a quasi- stationary state (QSS) in mechan- 
ical equilibrium. This is a stable stationary solution of the Vlasov equation (on the coarse- 
grained scale) that is usually not described by the Boltzmann distribution. On longer timescale 
tR ~ (N/ \nN)trj the encounters between stars (due to finite N effects) have the tendency to 
drive the system towards a statistical equilibrium state described by the Boltzmann distribution. 
In reality, this process is hampered by the escape of stars and the gravothermal catastrophe. 
The collisional evolution of the system is described by the Landau-Poisson system which is the 
1/N correction to the Vlasov limit (it singles out the Boltzmann distribution among all sta- 
tionary solutions of the Vlasov equation) [T2] . In fact, due to the time scale separation between 
the phase of violent relaxation (inertial effects) and the phase of collisional relaxation (finite 
N effects), we can consider for intermediate times that the distribution function is a quasi- 
stationary solution of the Vlasov equation of the form / = f(e,t) (for spherical systems) that 
slowly evolves under the action of close encounters according to the orbit-averaged-Landau 
equation (traditionally called orbit-averaged- Fokker-Planck equation). This implies that the 
lifetime of the QSS is long as it increases as a power of N. It slowly evolves under the effect of 
encounters which act as a perturbation of order 1/N with N ^> 1. Therefore, the system first 
reaches a state of mechanical equilibrium (through violent relaxation) then a state of thermal 
equilibrium (through stellar encounters). These different phases of the dynamical evolution of 
Hamiltonian stellar systems have been studied by astrophysicists for a long time pQ. 

(2) The case tn <C to <C £r (£ — > +00) corresponds to the overdamped limit of the 
Brownian model. The velocities first relax towards the Maxwellian distribution on a timescale 
tn ~ £ _1 (due to the thermal bath) and the density relaxes towards a state of mechanical 
equilibrium on a longer timescale (Smoluchowski diffusive time). Therefore, the system first 
reaches a state of thermal equilibrium (because of the terms of friction and noise in the Langevin 
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equations) then a state of mechanical equilibrium (through inertial effects). This overdamped 
regime, described by the Smoluchowski-Poisson system, has been studied in our series of papers 

E El DEI El HDI El- 

(3) Finally, there is an interesting case tp <C t B <C t R that has not yet been studied. In 
that case, there is first a phase of violent relaxation on a timescale ~ tp leading to a quasi- 
stationary state (QSS) in mechanical equilibrium like in case (1). This phase is followed by 
a thermalization leading to the Boltzmann distribution on a timescale ts ~ due to the 
thermal bath, i.e. the combined effect of imposed friction and diffusion in the Langevin model 
(HI)-©, n °t to "collisions" (finite N effects) as in case (1). The first phase is described by 
the Vlasov-Poisson system and the second phase by the Kramers-Poisson system. For £ — > 
(but £ > (lniV/iV)^ 1 ) there is a time scale separation between the phase of violent relaxation 
and the phase of Brownian relaxation. Similarly to case (1), we can consider for intermediate 
times that the distribution function is a quasi-stationary solution of the Vlasov equation of the 
form / = /(e, t) (for spherical systems) that slowly evolves under the action of imposed friction 
and diffusion (thermal bath, not collisions) according to the orbit-averaged-Kramers equation 
derived in Sec. 12.61 Since the Brownian timescale ts ~ is independent on N, this implies 
that the lifetime of the QSS in this regime is independent on N. Furthermore, it is shorter 
than in case (1) if £ 3> (In N/N)tp . Therefore, the system first reaches a state of mechanical 
equilibrium (through violent relaxation) then a state of thermal equilibrium (through the effect 
of imposed fluctuation and dissipation, i.e. the thermal bath). This is the opposite situation to 
case (2). The study of the orbit-averaged-Kramers equation will be considered in a future work. 
Note that if ts and £r are comparable, one must take into account simultaneously the effect of 
the thermal bath (friction and random force) and the effect of collisions (finite N effects). This 
is another interesting case. Finally, we stress that these different regimes should be observed 
for other potentials of interaction u(\r — r'|) than the gravitational one (e.g. for the HMF and 
BMF models j33j). Kinetic theories of Hamiltonian and Brownian particles with long-range 
interactions are discussed in [T2| at a general level. 

Acknowledgements We are grateful to the referees for useful remarks that helped to 
improve the presentation of the paper. 

A Eigenvalue equation for polytropes 

In the case of a polytropic equation of state, we can put the eigenvalue equation (pjj) in a 
dimensionless form |2~6] : 

d (e x ~ n dq 

di\^Ti 

where in the Emden function [34J and 

n 2 

is the dimensionless eigenvalue. This relation shows that the dimensional eigenvalue scales as 
A(A + £) oc Gp where p is the central density. Now, for a polytrope, the central density is 
connected to the mean density p by a proportionality constant depending only on the polytropic 
index n [3I|. Therefore, we also have A(A + £) oc Gp. 



+ 



nq nfl? 



A(A + Q 
SdGpo 



-1 



(95) 



(96) 
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B The exact Virial theorem in d dimensions 



We consider iV Brownian particles with mass m a in gravitational interaction. Their equations 
of motion are 

*? = £ G 7 {x * "if } - + y/™:m). (97) 

Here, the greak letters refer to the particles and the latin letters to the coordinates of space. 
For simplicity, we have assumed that the friction coefficient £ is the same for all the particles 
but this assumption can be relaxed easily (see below). The diffusion coefficient is given by the 
Einstein formula D a = £kBT/m a . The multi-species Smoluchowski-Poisson system has been 
studied in [11 . In this Appendix, we establish the exact Virial theorem associated with the 
stochastic equations ([9*7]). The moment of inertia tensor is defined by 

h, = J2 m - x ? x r ( 98 ) 

a 

We introduce the kinetic energy tensor 

K l3 = \Y. m « x * x l (") 

a 

and the potential energy tensor 



Wij = (?V m a mp 

1 (sf-sfXtf-sf) 

■-G^m a mp J — — J —, 100 

2 ^— ' r« — r„r 



where the second equality results from simple algebraic manipulations obtained by exchanging 
the dummy variables a and (3. Taking the second derivative of Eq. (|98jl. using the equation of 
motion (|97j). and averaging over the noise and on statistical realizations, we get 

7^4) + ^(4) = 2(A- y > + - i ^"(P^ + P^) dft, (101) 

where we have included the Virial of the pressure force Pj = Pn-AS^ on the boundary of a box: 

(4W = E 6o , EJ^ + ^X) = Z box Ea( p ?k x f + p j a k x ?) AS k 3 - The scalar Virial theorem 
is obtained by contracting the indices 

-(/) + -£(/) = 2(K) + (W u ) - jp ik XidS k , (102) 

where 

I = ^2m a r 2 a , K = ^^2m a v 2 a , (103) 



3 Note that if the particles have a different friction parameter £ a , the term £/ is replaced by ^2 a ^ a m a (ii xf + 
xfx") which can also be written ^ s where (s) denotes the different species of particles and Zy )S is the 

moment of inertia due to particles of species s. 
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are the moment of inertia and the kinetic energy. On the other hand, 



jo 1 P 1 n 1 



For d 7^ 2, we find that 



Wi* = (d - 2) W (105) 



where W is the potential energy 



W = m ° m * (loe) 

2(d-2)^|r,-r a |^ 1 Ubj 

In that case, the scalar Virial theorem reads 

^(7) + ^£(7) = 2(K) + (d - 2)(W) - jp ik XidS h . (107) 

For Hamiltonian systems {D — £ = 0), the total energy E = K + W is conserved. Hence, the 
Virial theorem can be written in an unbounded domain (P = 0): 

-I = 2K+(d-2)W = 2E + (d-A)W. (108) 

This is the extension of the usual Virial theorem in d dimensions (this equation is exact without 
averages). We note that the dimension d = 4 is critical as also noticed in [35] using different 
arguments. In that case, I = AE which yields after integration I = 2Et 2 + Cit + C2- For E > 0, 
I — > +oo indicating that the system evaporates. For E < 0, I goes to zero in a finite time, 
indicating that the system forms a Dirac peak in a finite time. More generally, for d > 4 since 
(d — 4)W < we have / < 2Et 2 + C\t + Ci so that the system forms a Dirac peak in a finite 
time if E < (this remains true for a box-confined system with a pressure force at the wall). 
Therefore, self-gravitating systems with E < are not stable in a space of dimension d > 4. 
The study in Ref. jHS] indicates that this observation remains true even if quantum effects 
(Pauli exclusion principle) are taken into account. This is a striking result because quantum 
mechanics stabilizes matter against gravitational collapse in d < 4 [36] . For 2 < d < 4, since 
(d — 4)W > we conclude, according to (|108j) . that the system evaporates if E > while an 
equilibrium is possible (but not compulsory) if E < 0. Finally, for d < 2, since W > 0, the 
energy is necessary positive (E > 0). Since (d — 4)W < an equilibrium state is possible. 
For d = 2, we have the simple result 



—Gj2 m *™p- ( 109 ) 



For equal mass particles, 



which reduces to 



W it = ~GN(N-l)m 2 (110) 



W„ = -^ (111) 
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for N ^> 1. We also note that 

Y,™ a mp = M 2 -J2 m l ( 112 ) 

The first term is of order N 2 m 2 and the second of order Nm 2 (where m is a typical mass). 
Therefore, in the mean-field limit N — > +00, we recover Eq. (jlllj) even for a multi-components 
system, cf [TT] . 

At equilibrium, the Virial theorem (jl(J2|) reduces to 

2{K) + {W u ) = j P ik XidS k . (113) 

If the system is at statistical equilibrium, then (K) = ^Nk B T and P^- = pdij with p = 
Yls Ps^BT/m s [TT]. Then, introducing the notation (26) of Paper I, we get 

dNk B T + (Wu) = dPV. (114) 

For an ideal gas without self-gravity (Wu = 0), we recover the perfect gas law PV = Nk B T. 
Alternatively, for a self-gravitating gas in two dimensions, using Eq. (jlllj) . we get the exact 
equation of state 

PV = Nk B (T-T c ), (115) 

with the exact critical temperature: 

AN K J 



For equal mass particles, we get 



Gm 2 

k B T c = —(N-l), (117) 



as in Eq. (j77j). In the mean-field limit 

GM 2 . s 

k B T c = . 118 

AN K J 

We now consider the strong friction limit £ — * +00. To leading order in l/£, the iV-body 
distribution is given by 

P N (r u v a , r N , v N , t) = e-^E^i ^^(r!, rjv , t) + 0(1/0, (119) 

as can be deduced from the multi-species iV-body Fokker-Planck equation ^1] generalizing Eq. 
P)l [this distribution is obtained by requiring that the Fokker-Planck collision term remains 
finite when D, £ — > +00]. From this expression, we find that P^ = with p = p s k B T/m s 
and (Kij) = ^Nk B T6ij even for the out-of-equilibrium problem. From Eq. (jlOlj) . we obtain 
the overdamped Virial theorem for a self-gravitating Brownian gas 

~£(4-) = (Wij) + Nk B TSij - ^j>p(8 ik x j + S jk Xi) dS k . (120) 

We can obtain this result from a different manner. In the strong friction limit £ — > +00, the 
inertial term in Eq. (j§Tj) can be neglected so that the stochastic equations of motion reduce to 

A? = W £ G 7 ( -r"|f ) + (121) 



lr/3 - r r 
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where D' a = k B Tp a and p a = l/£m a . Taking the derivative of the tensor of inertia (|9*HJl and 
using Eq. (jl21j) . we get 

4- = jWn + m a ^/2D^[xfR'»(t) + x<*R<*{t)\. (122) 

a 

Now, averaging over the noise using (x^R?) = ^2D' a 5ij, and on statistical realizations, we 
recover Eq. (|12U|) . The scalar Virial theorem reads 

-£(/) = dNk B T + (W u ) - dPV, (123) 

where P is defined as in Paper I. This is similar to Eq. (32) of Paper I obtained from the SP 
system but now Wu is given by Eq. (jl()4j) . In particular, in d = 2, using Eq. (jl()9|) we obtain 

^£(7) = 2Nk B (T - T c ) - 2PV, (124) 
with the exact critical temperature ()116j) . 



C Dynamical stability of homogeneous systems 

In this Appendix, we study the linear dynamical stability of a stationary solution of the damped 
barotropic Euler equations (|2^ j) - (j2^j) that is infinite and homogeneous, i.e. Po( r ) = P an d u = 0. 
For sake of generality, we consider a potential of the form 

$(r,t) = J u(r-r')p(r',t)dr', (125) 

where u(\r — r'|) is an arbitrary binary potential of interaction. We shall thus obtain a gener- 
alization of the Jeans instability criterion. We note that an infinite homogeneous medium is 
a stationary solution of the damped barotropic Euler equations provided that it satisfies the 
condition of hydrostatic balance Vpo + PoV$o — which reduces to V$o — 0. With Eq. (J125J) 
this can be written 

J ~Q~dx = or J u(x)dx. < oo. (126) 

We shall assume that this condition is fulfilled. We note that for the gravitational potential, 
this condition is not fulfilled since V • / Vwrfx = f Awc/x = S^G J 5(x)c?x = S^G ^ 0. Still, 
the equations for the perturbation are well-posed mathematically and, neglecting the above- 
mentioned inconsistency at zeroth-order, can be considered as a first step to investigate the 
dynamical stability of a gravitational system. This is the so-called Jeans swindle pQ. The 
linearized damped barotropic Euler equations can be written 

^ + pV-5u = 0, (127) 



p— = -cjA5p - pA5$ - £p5u, 



'1281 



5$(r,t) = / u(t -r')5p{r',t)dr' 



(129) 
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where we have introduced the velocity of sound c 2 = p'(p). They can be combined to give 

^£ + ^ = 4A Sp + P A6Z, (130) 
with 5<£> = u * 5p. Looking for solutions of the form 

5p(r,t) = J 8p{k,u)e i( - 1 '- r -" t) dkduj, (131) 

we obtain the general dispersion relation 

u(uj + it) = c 2 s k 2 + (2n) d u(k)pk 2 . (132) 

Introducing A = — iuj and v(k) = —(27r) d u(k), this can be rewritten 

A 2 + £A- k 2 {pv{k) -c 2 ) = 0. (133) 

The solutions are A± = §(-£ ± VA) with A(k) = £ 2 + Ak 2 (pv(k) - c 2 ). If c 2 < pv(k), then 
A(A;) > £ > and the system is unstable as A + = §(— £ + v^A) > 0. If c 2 < pv{k), either 
A(k) < and R e (X) = — ^/2 or < A(k) < £ implying A± < 0, so the system is stable. 
Therefore, the system is stable if 

c 2 s > pv(k), (134) 

and unstable otherwise. For attractive potentials v(k) > 0, this gives rise to the existence of a 
critical point as discussed in [TJ|. Indeed, a necessary condition of instability is that 



pv{k) max . (135) 

If this condition is fulfilled the range of unstable wavelengths is determined by 

v{k) > c 2 Jp, (136) 

and their growth rate is \+(k). For the Euler equations (£ = 0), the dispersion relation reduces 
to 

uj 2 = c 2 s k 2 + (2n) d u(k)pk 2 . (137) 

In the unstable case the perturbation grows exponentially and in the stable case the per- 
turbation presents undamped oscillations. In the overdamped limit of the model £ — > +oo 
(Smoluchowski), the dispersion relation reduces to 

i& = c 2 k 2 + {27i) d u{k)pk 2 , (138) 

which can be obtained directly from the generalized Smoluchowski equation (|35j). In the unsta- 
ble case the perturbation grows exponentially and in the stable case the perturbation decreases 
exponentially. The intermediate case of finite friction can be treated as in Sec. E]but explicit 
results demand to specify the potential of interaction. 

For the attractive Yukawa potential v(k) = SdG/{k 2 + k 2 ) jT2j, we find that the system is 
unstable if c 2 < (c 2 ) cr i t = SdGp/k 2 . for the wavevectors such that 

k < k max = \ 5 k 2 ,. (139) 

V C s 
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The growth rate \+(k) is maximum for k 2 = (SdGpk 2 /c 2 ) 1 / 2 — k\ and its value A* = A + (&*) is 
given by 



2A* = -Z+J? + AS d Gp (l - ^c 2 k 2 /S d Gp^ . 



(140) 



For c s = (cold systems), we have k max = k* = +00 and A* = ~(— £ + a/£ 2 + ASdGp). 
The system is most unstable at small wavelengths. For k = (gravitational potential), we 
have {c 2 s ) crit = +00, k max = kj = (S d Gp/c 2 s ) 1/2 (Jeans length), k* = and A* = + 
a/£ 2 + ASdGp). The system is most unstable at large wavelengths. For £ = (Euler), 
A* = {SdGp) 1 ' 2 ^ - ^/c 2 k 2 /S d Gp) and for f -> +00 (Smoluchowski), A, = (SdPp/0(l ~ 
a/ c 2 s k^/ SdGp) 2 . Other examples of potentials of interaction are given in 0121 . 
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